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Abstract 

We consider interacting Bose gas in thermal equilibrium assuming 
a positive and bounded pair potential V{r) such that < / drV{r) = 
a < oo. Expressing the partition function by the Feynman-Kac func- 
tional integral yields a classical-like polymer representation of the 
quantum gas. With Mayer graph summation techniques, we demon- 
strate the existence of a self-consistent relation = F{fj. — ap{fj,)) 
between the density p and the chemical potential fi, valid in the 
range of convergence of Mayer series. The function F is equal to 
the sum of all rooted multiply connected graphs. Using Kac's scaling 
V^(r) = 7^1/(7r) we prove that in the mean-field limit 7—5-0 only tree 
diagrams contribute and function F reduces to the free gas density. 

We also investigate how to extend the validity of the self-consistent 
relation beyond the convergence radius of Mayer series (vicinity of 
Bose-Einstein condensation) and study dominant corrections to mean 
field. At lowest order, the form of function F is shown to depend 
on single polymer partition function for which we derive lower and 
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upper bounds and on the resummation of ring diagrams which can be 
analytically performed. 
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1 Introduction 

The interest for a better understanding of Bose-Einstein condensation has 
been strongly stimulated in recent years by the beautiful experimental ob- 
servations of condensates of cold atoms in traps [1,2]. 

Concerning rigorous results on the existence of Bose-Einstein condensa- 
tion in an interacting gas with pair interactions, we quote the work of Licb 
and Seiringer [3] . The authors show the existence of off-diagonal long-range 
order in the ground state of a system of Bose particles confined by an external 
potential in the dilute limit in which the Gross-Pitaewski equation becomes 
exact. Siito shows condensation for the trapped gas at non zero temperature 
[4]. In [5] Lauers, Verbeure and Zagrcbnov prove the existence of a Bose- 
Einstein phase transition for the homogeneous gas under the assumption that 
there is an energy gap at the bottom of the one-particle spectrum. However, 
to our knowledge, there is still no proof of Bose-Einstein condensation in 
the interacting gas when there is no trap and no gap. In the present work 
we revisit this venerable many-body problem with a new point of view, the 
technique of quantum Mayer graphs. 

In order to provide orientation and motivation for our approach we recall 
some facts pertaining to the mean-field ^ Bose gas at an heuristic level. The 
Haiiiiltoiiiaii of Bose ])articles in a volume A interacting with a constant 

""^Here we use "mean field" in the sense of van der Waals, that is collective effects of 
a long range interaction. The same word is also used in the context of the dilute limit 
[3] leading to the Gross-Pitaewski regime. In fact, in the latter case, the situation is the 
opposite: dominant effects of the interaction are due to rarefied local binary collisions. 
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repulsive potential of strength a/|A|, a > 0, reads 



R + ^ N{N-l) 

J^mi,N — J^O,N + 1^ [i-) 

where Hq n is the total kinetic energy. It leads to the free energy density 

/mf(/?,p) = /o(Ap) + |p' (2) 

where p is the particle density, /3 is related to the temperature T by /3 = 
1/kBT, ks denotes the Boltzmann constant, and fo{P,p) is the free energy 
of the non-interacting Bose gas. Differenciating with respect to p yields the 
relation 

/^o(^, p) = Pmf p) - ap (3) 

between the chemical potential po{P,p) of the free gas and the chemical 
potential /imf(/?,p) of the mean-field gas as functions of the density. Since 
the grand canonical densities po of the free gas and the density pmf of the 
mean-field gas, considered as functions of the respective chemical potentials, 
are the inverse functions of po{P, p) and PmiiP, p) at fixed /3, Q is equivalent 
to the self-consistent equation 

Pmf p) = Po(/5, P - apmf p)) (4) 

In (gj) 

^o(^' ^) = (27rA2)3/2 ^ 

which is the well known formula for the grand canonical density of the free 
gas with 



A = h^^J^ (6) 



representing the thermal wavelength [6]. 

The series (jS)) converges for /i < 0, so that the self-consistent equation 
dH) is meaningful whenever u = p — apmf p) < 0, namely for p < pc where 
the critical chemical potential pc is given by 

Pc - apnii{(3,pc) = 0. 

Thus Pc has the value pc = apo,c(/5) where Po,c(/5) = Po(/^)P = 0) is the 
critical density of the free gas. 
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At this point it is worth noting that the solution Pmf(/^) of (jH) for /i < /ic 
can be extended to the range > /Xc by continuity. Indeed, we differentiate 
equation Q with respect to the chemical potential to obtain ^ 

One sees from (jSJ that {pq)'{u) = oo for z/ > 0. Hence (pmf)'(/^) = 1/'^; z/ > 
and requiring the continuity of the density at /i = pc gives 

Pmf(p) = ^, P > Pc (8) 

Equations (jH) and (jHl) define the density for all values of the chemical poten- 
tial. For ji > fic, there is a Bose condensate of density Pmf(p) — Po,c- These 
facts have been established with full mathematical rigor in several works [7, 
8, 9] and see [10] for a review. In particular they are obtained for the Kac 
interparticle potential 

V,{t) = 7V(7r) (9) 

in the scaling limit of 7 tending to zero. When 7^0, the potential Vy{r) 
extends its range to infinity whereas its amplitude tends to zero in such a 
way, that the mean potential energy 

J dr-f^VM = a (10) 

stays constant. It is well known that in classical statistical mechanics this 
is the appropriate limit to rigorously recover the van der Waals mean field 
theory [11]. Also the methods of Mayer graphs enable to calculate corrections 
to the mean field limit for small 7 [12]. 

In this paper we propose a similar approach to study an interacting Bose 
gas with non-singular repulsive interactions. In section 2 we recall the "poly- 
mer" representation of the Bose gas in thermal equilibrium. Combining the 
Feynman-Kac functional integral representation of the Gibbs weight together 
with the decomposition of permutations into cycles, one finds the grand- 
canonical partition function in the classical-like form of a gas of interacting 
polymers. Polymers are Brownian closed loops associated with a number of 
Bose particles belonging to a permutation cycle. Each loop has a self en- 
ergy and there is a loop-loop pair potential. In the space of polymers, all 

^We keep the temperature fixed and omit from now on (3 in the notation. 
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the techniques of classical statistical mechanics are available, in particular 
the analysis of the partition function and of the density with the Mayer dia- 
grammatic techniques. In this way we show in section 3 how the mean field 
equation Q can be recovered by summing the tree graphs. 

In section 4 we treat an interacting gas with a general short range repul- 
sive two-body potential and establish that its density obeys an exact equation 
of the form 

p{fi) = F{fi - ap{fi)) (11) 

The function is defined as the sum of multiply connected Mayer graphs. 
We prove with the help of the Penrose tree-graph inequality that the corre- 
sponding diagrammatic expansion of F is convergent at low density, namely 
if the chemical potential is sufficently negative. Equation (jllj) provides a 
generalization of the mean-field equation to the interacting gas: indeed 
if one introduces a Kac potential V^(r), it is seen that F^{fj,) reduces to the 
density po{fi) of the free gas as 7 ^ 0, thus recovering equation (jD). The 
central goal then is to extend the validity of equation ()11|) to higher densities 
(hopefully up to a critical density), that is, to approach the Bose transition 
point from the dilute phase, as in the mean field theory. 

In section 5 we discuss the mathematical problems that arise at this point. 
The study of the critical point will require the control of the asymptotic be- 
haviour of the partition function of a single long repulsive polymer as well as 
the mutual interactions between different polymers. One should note how- 
ever that the polymers occuring in the representation of the Bose gas differ 
from the standard self-repelling classical polymers because of the specifically 
quantum mechanical "equal time interaction" introduced by the Feynman- 
Kac formula. Hence, the results of the theory of classical polymers cannot 
be used without further consideration. As first investigation in this direc- 
tion, we give lower and upper bounds on the partition function of a single 
polymer indicating that equation (fTT|) continues to hold for a range of chem- 
ical potentials larger than that assuring the convergence of the Mayer series. 
Moreover we can calculate the effects of mutually interacting polymers at 
lowest order in the Kac parameter 7 by summing up the ring diagrams in a 
closed form. Concluding remarks are presented in section 6. Proofs of some 
lemmas are relegated to appendices 1 and 2, and appendix 3 is devoted to 
the extension of our methods to an inhomogeneous Bose gas confined by an 
external potential. 
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2 The polymer representation of the Bose gas 



We consider bosons of mass m in three dimensions with Hamiltonian 

HN = -^j:A,+ E nr.-r,) (12) 

enclosed in a box A, A being the Laplacian with Dirichlet conditions at the 
boundary of A. The pair potential V{r) is short range, repulsive and without 
singularities: 

V{r) > 0, J drV{r) = a < oo, sup V{r) = V <oo (13) 

We assume moreover that its Fourier transform V^(k) is positive. 

The polymer representation of the grand-partition function at chemi- 
cal potential /i and inverse temperature /3 arises when the Gibbs statisti- 
cal weight is expressed in terms of the Feynman-Kac path integral where 
quantum fluctuations are represented by Brownian trajectories. The open 
trajectories associated with exchange contributions are reorganized in larger 
closed loops (or polymers) containing several particles according to the de- 
composition of permutations into cycles. The result is that the grand parti- 
tion function of the Bose gas can be written in a classical-like form as (the 
so-called magic formula) 

oo 1 „ n 

= E - / n dC.ziC,) exp[-/?[/(£i, . . . , £„)] (14) 

provided that suitable definitions of the phase space integration and of the 
interaction are given. In one form or another, this representation has been 
known since a long time in various contexts starting with the work of Ginibre 
on the convergence of quantum virial expansions [13]. It is also used to 
implement numerical simulations of the Bose gas [14]. The present form 
(HD) has been derived and applied by Cornu [15] to Coulomb systems and we 
follow here the definitions given in Chap.V of [16]. A self contained derivation 
can also be found in [17]. 

An element £ of the phase space, called a loop or a polymer, 

C = (R, q, X(s), < s < g) (15) 
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Figure 1: A g particle loop 



is specified by its position R in space, the number q of particles belonging to 
it and its shape X(s). The q particles are located at positions 

Tfc = R + AX(A;- 1), A; = l,...,g, q + l = l (16) 

and 

rfc,fe+i(s) = R + AX(A;- 1 + 0<s<l (17) 

is an open path joining the k particle at to the k + 1 particle at r^+i where 
A is the thermal wave length (jH)). The loop can be viewed as an extended 
object at R that has internal degrees of freedom (g, X) with q the number of 
particles belonging to a permutation cycle and X the shape of the loop, see 
Fig. 1. 

The thermal wavelength A (the only place where the Planck constant 
occurs) gives the extent of quantum fluctuations. The shape of the loop X is 
a Brownian bridge (a closed Brownian path), parametrized by the "time" s 
running in the "time interval" [0, g] with X(0) = X(g) = 0. It is distributed 
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according to the normalized Gaussian measure / Dg(X.) . . ., with covariance 



J D,{X)X,{s,)X,{s,) = 6,,q 









mm — , 


-) 








q q 



(18) 



where = 1,2,3 are the Cartesian coordinates of X. Integration on 

phase space means integration over space and summation over all internal 
degrees of freedom of the loop 

J dC--- = f^J Dg{X.) J^dR--- (19) 

Because of the Dirichlet boundary condition, the paths X are constrained to 
stay in the volume A, but we do not write this constraint explicitly since it 
will be removed later in the infinite volume limit. 

The interaction energy of two loops Ci, Cj is the sum of pair Feynmam- 
Kac potentials between the particles associated to the loops 



V{Ci,Cj) = f dsV(R^ + \Xi{k-l + s)-Rj-\Xj{i-l + s)) 

k=l £=1 

= / dsi dsj6{si - Sj)V{Ki + AXj(si) - R,- - XXj{sj)) 
Jo Jo 

(20) 

In fpUj) . the distribution 

oo 

S(s) = e^'™ (21) 

n=— oo 

is the periodic Dirac function of period 1. 

The activity z{C) of a loop is related to the chemical potential /i of the 
particles by 

'^^^ = q{2nqX^r/^ exp(-/3^(/:)), z = e^^ (22) 

It incorporates the interactions U{C) of the particles in the same loop (the 
self energy of the loop) 

U{C) = i r dsi r ds26{si - S2)V{X{X{si) - X{S2))) - ]-qV{0) 
Z Jo Jo 2 

(23) 
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The last term substracts out the self energy of the particles. U{C) can as 
well be written as 



f/(^) = 2^'^^i/J^^25(3i-S2)(l-5f,'i,[,,])nA(X(s0- >0 (24) 

making manifest that U{C) is a positive quantity (in the Kronecker delta 
symbol [s] denotes the integer part of s). As a direct consequence one gets 
the bound 

From the structure ()14|) of the partition function and the above defini- 
tions, it is clear that the calculation rules of classical statistical mechanics 
apply to the system of loops. We shall take advantage of this fact to analyze 
the system of loops first and derive from there the results for the original 
quantum gas of particles. In particular all the powerful techniques of Mayer 
graphs are available to expand the loop density and the loop correlations in 
powers of the loop activities z{C). It is convenient to introduce the abbrevi- 
ated notation Ci = i and dCi = di. Mayer bonds are defined by 

/(z,j) = e-'^^(^'^)-l (26) 

and weights at vertices by z{i) ((221) • Integration at vertices di has to be 
performed according to (fT^ . 

Notice that the bond is integrable over space since from the positivity of 
y(r), 

|/(2,j)|< (27) 

and from ((^0]) 

JdR,\f{C„q)\ < (3jdR,V{C,,C,) 

= (3 P dsi p dsj5{si - Sj) J dRV(R) 

= f^aqiqj (28) 

The loop density pioop('C) is given by the standard expansion 

oo 

Ploop(l) = E^^n(l) (29) 

n=l 
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where 



BJl) = y — d2--- dnz(l)z(2) ■ ■ ■ z(n)u(l, 2 




(30) 



n=l 



and 



«(l,2,...,n)=^ n 



fi^,J) 



(31) 



is the Ursell function. The sum runs over all labelled connected graphs r„ 
with n vertices. We have directly written the loop density (jHUj) in the infinite 
volume limit, namely extending the spatial integrals (iR2, . . . (iR„ at vertices 
2, . . . ,n over the whole space, whereas the vertex 1 (the root point of the 
graph) carries no spatial integration. The existence of the infinite volume 
limit of individual Mayer graphs follows from the translation invariance and 
integrability of the Mayer bond as in the classical case. As a consequence 
Pioop('^) = Pioopil,'^) does not depend of the location R of the loop C in 
space. 

Finally, to obtain the original particle density p(/i) from the loop density, 
we have to sum pioop(jC) over the internal variables of £, 



the aditionnal q factor taking into account that the loop C carries q particles. 

Although classical methods have been used, it is important to stress the 
difference between the loop representation of the equilibrium state of the 
quantum Bose gas and that of a gas of genuine classical polymers. First the 
chemical potential /i is not the variable thermodynamically conjugate to the 
polymer number, but to the original particle number. Moreover, the loop 
interactions differ from the classical polymer intractions by the quantum 
mechanical "equal time prescription" which originates from the Feynman- 
Kac formula. This equal time prescription is manifested in (j^UI) and ()23p 
by the occurence of the periodic delta function 6{si — S2). In the classical 
interaction, every segment of a polymer interacts pairwise with any other 
segment, which would correspond to the potentials (j20p and (j23|) without the 
equal time prescription. The purely classical gas of point particles interacting 
by means of the two-body potential V^(r) is recovered if the thermal lenth A 
is set equal to zero and if only terms with q = I (Boltzmann statistics) are 
retained. 




q=l 



(32) 
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Figure 2: The bonds f^^^i,]) and f^^\i,j) 

3 Mean field limit and tree graph summation 

Let us first show how the self-consistent mean field Bose gas is recovered in 
the diagrammatic analysis. To this end we split the bond f{i,j) into its part 
linearized in the potential and higher order terms 

f^'\^,J) = e-^^(^'^)-l + W^,j) (33) 
These bonds are represented in Fig. 2. 



This simply enlarges the previous class of graphs F to graphs (still denoted 
by F) where each bond can be either f^"'\i,j) or f^''\i,j). 

In the mean field limit 7 — > 0, the contributions of order 0{1) in 7 due 
to the scaled potential 7^\^(7r) will come from the linearized bonds 

4'^)(2,j) = -13^^ £ds, |J'rfs,5(s,-s,)V(7(R.-R,)+7A(Xi(s,,)-X,(s,))) 

(34) 

and from the vertex weights 

The bond f!f'\i,j) will give contributions of order 0{1) since under scaling 
fllO|) its total spatial integral 

J dRjf^''\C,,,) = -I3aq,q, (36) 

is independant of 7. In the activity (f^^ we simply disregard the self-energy 
U^{C) since the latter is 0(7^). 
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''"n+1 



Figure 3: Rooted tree T„_,.i with subtree T, 



Proposition 1 

The density ptree(/^) calculated as the the sum of all tree graphs with bonds 
and vertices U^) verifies the mean field equation 

Consider a rooted tree graph T„_|_i with vertices (0, 1, . . . , n) for which the 
root point is of degree 1 (0 is the label of the root point) ^, Fig. 3. Hence the 
root point is linked to the rest of the graph by a single bond, say /'•'^^(0, 1). 
Call T„ the subgraph of Tn+i with vertices (1,2,..., n) and t„(l) the value 
of this subgraph once integrated on the vertices 2, . . . , n. Then the value 
tn+i(>Co) of the rooted graph T^+i is 

t„+i(/:o) = -/5^^°H?o) / dC^V,{Co,CMl^i) (37) 



The root point can be linked to any of the n vertices of the graph T^. 
The resulting tree graphs are different beacuse they contain different links 
with the root point. Moreover, because of translational invariance tn(l) = 
t„(gi,Xi) does not depend on the position Ri of the loop Ci. Thus the 
spatial integration over Ri can be performed on V^{Cq, Ci) as in so that 
the total contribution reads 



oo „ 

Y^qJ D(Xi)nt„(gi,Xi) 

qi = l 



^The degree of a point is the number of lines incident at this point. 



(38) 
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According to equations and (jH^ the quantity in the bracket is precisely 
the contribution to the particle density ptree(/^) of the graph T„. Therefore the 
sum of all tree graphs having a root point of degree 1 is z^^^ (go) [— /5agoPtree(/^))- 
The sum of all tree graphs with root point of degree n is z^^^ (go) (— /5agoPtree(/^)]"'/ n\ 
(the factor l/n\ takes care of the fact that the labelling of vertices belonging 
to different branches attached to the root point can be permuted without 
giving rise to new Mayer graphs). Finally, summing on all trees rooted at 
the point gives the density of loops. According to (1221), obtain the par- 
ticle density we still have to sum on the internal variables go, Xo of the root 
loop with a factor go. Hence using we find 

Ptree(/^) = £ X] ( "/^"^OPtree (/^) (" 

go=ln=0 ^■ 

_ exp(/3(/i - aptrec(/i))) 

- ^j^, (27rgoA2)3/2 

= po(p - aptrcc(p)) (39) 
which is the mean field equation (^j). 

We now come back to the general Mayer series ()3U|) and prove its conver- 
gence at low density. 

Proposition 2 

The Mayer series / T^) converges for fi < — apo,c with po,c = Po{P,fJ^ = 0) the 
critical density of the free gas. 

We use the Penrose tree graph inequality for positive potentials which 
states that the the sum of the tree graphs provides an upper bound for the 
Ursell function [18] 

u{l,2,...,n) < E n 

T„ (i.j)GT„ 

< E n im^j)) (40) 

Tn (j.j)er„ 

the second inequality being a consequence of (j2Zj)- Moreover, if in the series 
finn|) we use the inequality we see that the tree summation with bonds 
|/(")(i,j)| = pv{i,j) and vertices z^'^\i) provides also an upper bound. Ac- 
cording to the Proposition 1 the latter series sums up to a function p(/i) that 
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obeys 



p{fi) = po(/i + ap(/i)) 



(41) 



This is the mean field equation of a Bose system for the negative potential 
—V{r): It has a finite solution provided p + vp{fi) < which is equivalent 
to fi < — /?apo,c- We conclude that 

p(/u) < p(/u)) < oo, p < -apo,c (42) 

Notice that in the case of the scaled potential Vy{r), the convergence of the 
series (IHU)) defining p-y(p) is not only absolute but also uniform with respect 
to 7 . This follows from the fact that in the evaluation of the tree-graph 
contributions one encounters only the integrated bonds 

J rfR,|/f (A,,)| =/5ag,g, 

The Penrose inequality yields thus a 7-independent upper bound. 

Now we show that in Kac limit, the density converges to the mean-field 
value. 

Proposition 3 

Let p^(p) be the density associated with the scaled potential V^(r) = 7^V^(7r). 
Then 

linip^(p) = Pnif(/i), P < -apo,c (43) 
7— >o 

Consider first a graph r„ with vertices l,2,...,n (rooted at point 1) 
which is not a tree, i.e it contains at least one cycle. Delete some bonds in 
such a way that the graph thus obtained is a connected tree T„ with value 
tn{l)- In order to find an upper bound for the value gn,'y{^) of Tn we use the 
inequahtes 

\f,{i,j)\<PV,{t,j) (44) 
for the bonds remaining in the tree and 

\f,{t,j)\<l'PmV (45) 

for each deleted bond. Notice that (j45|) is obtained by replacing the pair 
potential in (jMj) by its supremum V. Moreover the activity z{i) is bounded 
by z^'^\i) (|^. Spatial integration of a bond {i,j) of the tree yields the factor 
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Paqtqj as before so that with a vertex (i) of degree fcj in r„ receives a 
factor gf' multiphed by the activity z^^\qi). If £ bonds have been deleted, 
the resulting inequality has the form 

oo 

\9nM\ < ii'vYiPar-' E Qi'---Qt^^'\qi)---^^'\qn) m 

q2,---,qn = l 

In view of the exponential factor e^^'^ in (j25p . the g-series in (I46|) is convergent 
for /i < 0. Hence the presence of one cycle in r„ implies that 5'n,7(l) tends 
to zero not slower than 7^ for 7 — > (and as 7^^ if r„ has i cycles). 
If r„ is a tree we decompose the Mayer bonds as in (jHHI) . Since 

J)l < I {V{^,J)C < l7'f3mV{(3V{t,j)) (47) 

all the contributions to r„ containing f^^\i,j) bonds vanish as 7 — 0. 

As we have seen from the Penrose estimate and from Proposition 2, the 
Mayer series constituting p^{fi) is absolutly convergent uniformly with re- 
spect to 7. By dominated convergence, the limit of the diagrammatic sum 
can be calculated term by term. Therefore we are left with the sum of trees 
where all the bonds are of the type f!f'\i,j) and where the activites z^{i) 
can be replaced by z^'^\i) as 7 — 0. This is exactly the situation of the 
Proposition 1 thus proving Proposition 3. 

4 The self-consistent equation for the inter- 
acting gas 

The reasoning which led to equation ()39|) was restricted to the tree diagrams. 
We now show how it can be generalized to provide an implicit equation for 
the exact density p(/u). To this end we consider the complete set of expanded 
Mayer graphs with the two types of bonds /*■"•'(«, j) and f^''\i,j) defined in 
and call f^^^^i ,j) a single interaction bond. A graph is said multiply 
connected ^ if it cannot be disconnected by cutting a single interaction bond. 
Then we define I{C) to be the value of the sum of all multiply connected 
graphs with one root point C, see Fig. 4. 



^In the context of Feynman diagrams, such graphs are also called one-line irreducibles. 
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i=o + 




+ 



Figure 4: The sum of multiply connected graphs 



Because of translation invariance, does not depend on the position 
R of £. Denoting here the chemical potential by , we define the function 



For V < —apo^c, the function F{i>) is represented by a convergent sum since 
it is a subseries of the absolutly convergent Mayer expansion (Proposition 2). 

Proposition 4 For p < —apo^c, p{fJ') verifies the equation 



To establish ()48j) and ()49j) . we proceed as in Proposition 1. Considering 
now a general structure of Mayer graphs we say that a vertex (i) in a graph 
has a star structure if there is a number of incident single interaction lines at 
(i) such that cutting any one of them disconnects the graph. Consider first 
a vertex (i) having a star structure consisting of one single interaction line. 
This interaction line, say /'•"^(«, 1), links the vertex (i) to a subgraph r„ with 
vertices (1,2, . . . ,n)] r„ has no other links with the rest of the graph. Fig. 
5. Once integrated on the points (2, . . . ,n), its value 5'n(l) does not depend 
on the position Ri of the loop Ci. Thus the spatial integration over Ri 
can be performed on the bond [—j3V{Ci, Ci)] yielding the factor {—(3aqiqi). 
Taking finally into account that the graph r„ can be attached through a 
single interaction line to the vertex {i) by any of its n vertices we find (as in 
dnZI)) that its total contribution to the vertex (i) equals 




(48) 



p(/i) = F{p - ap{p)) 



(49) 



z{Ci){-(3aqi) ^(li j F>q^{Xi)ngn{qi,Xi) 

qi = l 



(50) 
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Figure 5: Star structure at vertex (i) and subgraph r„ 



But from (jHzj) . the quantity in the square bracket is the contribution of 
the graph r„ to the exact density. Thus the sum of all such graphs con- 
tributes to the vertex (z) as z{Ci){—(3aqip{^)). Continuing the reasoning 
along the lines of the proof of Proposition 1 we consider next the sum all star 
structures at {i) having n single interaction lines. Their contribution equals 
z{Ci)[—l3aqip{ii)Y'/n\. Finally, summing all possible star structures at (i) we 
arrive at the formula for an effective activity 



z*{Ci) = z{Ci)e 



-f3ap(p)q^ 



g,(27rg,A2)3/2* 



(51) 



This is precisely the previous activity z{Ci) evaluated at the shifted 
chemical potential v = p — ap{p). 



lYiom. this analysis we see that p(/i) can as well be obtained by sum- 
ming all multiply connected graphs with bonds f^"'\i,j) and f^''\i,j) with 
effective activities z*{Ci) at vertices. Indeed, by construction, any graph in 
these series is an original Mayer graph. Conversely, each Mayer graph does 
appear therein. This can be shown by the following reasoning. Consider 
any connected graph made of the bonds f^"'\i^j) and f^^\i-ij) and remove 
all single interaction lines /"'{ijj) whose presence makes it not mutliply con- 
nected (by definition removal of such a line disconnects the root point from 
a part of the graph). Then the remaining subgraph connected to the root 
point is a multiply connected graph from which the original Mayer graph 
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is obtained in a unique way by forming star structures. This concludes the 
proof since, in view of formula (|51|). by forming F{i>) according to these rules 
we find the density equal to the value of this function at the shifted argument 
u = fi — ap{fi). 

As a corollary of Propositon 3, we have that if F^[v) is calculated with 
the scaled potential ^-(r) then lim^^o -^7('^) = Poi.^)- 

We can write I{C) appearing in the definition (PH|) of function -F^(z/) as 
the sum 

/(/:) = ^(/:)[i + /cycic(/:)] (52) 

of the pure root point term plus all the multiply connected graphs containing 
cycles and bonds. It is clear from the proof of Proposition 3 that 

the mean field density arises solely from the root point z{C) whereas /cycle (jC) 
becomes vanishingly small in the mean field limit. 

All the results presented so far are valid within the convergence radius of 
the Mayer series determined by < — apo,c- The question then arises if the 
density, as solution of (pUj) . can be extended to larger values of /z, as it was 
possible in the strict mean field case. An investigation of this question and 
of the mathematical difficulties involved is presented in the next section. 



5 The Bose gas beyond mean field 

5.1 Vertex contribution: the single polymer partition 
function 

In this section we study the structure of function F^p) in a more detailed 
way. The first term f^^^u) of the expansion of F{u) in multiply connected 
graphs corresponds to the root point z{C) which is a single vertex contribu- 
tion 

oo „ 
q=l 

= J^^2T.^<1) (53) 

Here we have introduced the mean value of the Boltzmann factor 

K{q) =< e~P^ J D,(X)e-^^('^'^) (54) 
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of a single closed polymer carrying the repulsive energy (j23|) . The Brownian 
bridge measure defining the average < ... >q is normalized to 1, namely if 
f/(g, X) = 0, one has n{q) = 1 and f''^\v) reduces then to the free gas 
density (0). Moreover, for positive potentials we have < ^(g) < 1. Clearly, 
the radius of convergence of the series ( ESI) is determined by the asymptotic 
behaviour of ^(g) for g ^ oo. The following bounds can be established 

Lemma 1 

There exists z/+ such that 

K{q) > 6"^''^+ (55) 
This lower bound follows from Jensen inequality 

K{q) > exp(-/5(f/),) (56) 

where 

{U), = jD,{^)U{q,^) (57) 

is the average energy of the polymer. The value of z/+ results from a direct 
calculation of {U)g presented in Appendix 1. 

There is also an upper bound: 

Lemma 2 

Let r be a fixed integer r > 2. There exists (depending on r) such that 
< < Z/+ and 

K[qr) < r'^/^ exp(— /3grz/_) (58) 

The upper bound is derived by splitting the closed loop X(s), < s < 
g, X(0) = X(g) = 0, into the union of two open Brownian paths Xi(s), < 
s < qi and X2(s), gi < s < g, and disregarding the (positive) interactions 
between the two paths Xi and X2. The details and the form of z/_ are given 
in Appendix 2. 

The lemmas imply that f^^\i^) is finite if z/ < z/_ and diverges for z/ > z/+. 
Let us suppose for a moment that there exists a critical value z/^ such that 

/(°)(^) < 00, z/<z/, 

/(°)(z/) = 00, u>u, (59) 

In this case, there will be a critical density pc = f^^\vc) and a critical chem- 
ical potential p^ = o-Pc + ^c- From the lemmas 1 and 2, one has necessarily 

z/_ < z^c < Z/_|_. 
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The determination of the possible existence of requires the knowledge 
of the exact asymptotic behaviour of K{q) as g — oo. In this respect let 
us make the following comment concerning the theory of classical polymers, 
defined as those interacting via the standard pairwise repulsion 

Ucx{C) = \r dsi r ds2V{X{sr) - X(S2)) (60) 

It is firmly established, although not rigorously proved, that the normal- 
ized partition function of a single classical closed polymer behaves (in three 
dimensions) as [19] 

nM = J /^,(X)e-^^-(^) ~ g ^3(.,.^i/2) ' 5 - ^ (61) 

where A is a constant depending of the choice of the potential V^(r) and 
t'poi = 0, 589 is the universal critical exponent for a swollen polymer. Inserting 
this asymptotic behaviour in the series gives Uc = A and a finite critical 
density pc since the series is convergent at = z/^. It is an open question 
to find out whether a similar situation holds for the "quantum" polymers 
subjected to the "equal time" interaction 



5.2 Bond contributions: interacting polymers 

Interaction between different polymers occurs in multiply connected graphs 
having bonds f^"'\i,j) or f^''\i,j)- It turns out that the subseries /ring('C) of 
-^cycie('^) defined as the sum of all multiply connected graphs having exactly 
one cycle of interaction bonds can be summed in a closed form, see Fig. 6 
(the first term of the series corresponds to the quadratic term V'^{i,j)/2 in 
the expansion of f^^\i,j), (jSSI))- The result is 



2 J h Jo 1 



-e 



2inn{s—t) 

(62) 

Here V^(k) is the Fourier transform of the potential and the positive coeffi- 
cients a^(k) 



k) = f^q Tds f Z},(X)z(/:)e'^'^-^(^)e2*™ (63) 

q=0 ''^ 



come from the summation on the internal degrees of freedom of loops at 
vertices. A derivation of formula (j62j) can be found in [20] where the effective 



20 



I = o + cOi + 

ring ^-^ ^-^*s,_^,^'^^ 



Figure 6: The sum of rings 



loop-loop potential is calculated as the sum of all chain graphs. Irmgi'C,) is 
given precisely by the formula (84) in [20] restricted to a single species of 
bosonic particles and with the Fourier transform 47r/|kp of the Coulomb 
potential replaced by the present short range potential V^(k) (in [20] the 
coefficients 47r/9Q;^(k) are noted K^(k,n)). 




We now check that /ring(>C) is also well defined in the range z/ < z/_ 
(or possibly u < ii there is a critical value Uc (j59|l ). As the modulus of 
the phase factors in (jH^ equals 1 the positivity of pV(k)a'^(k) implies the 
inequality 

^Hng(>C) < ^ / dk{(3Vm Y: «^(k) (64) 

^ n=-oo 

/,From dnSl) and ^ 

< f^'\i^) (65) 
We have defined 

oo „ oo pPvq 

/«(.) = ^«(^)^(^) = ^2.qX^Y/2 <<l) (66) 

and the inequality follows from the fact that the bracket [■ ■ ■] is less or equal 
to q . Hence the contribution of the ring diagrams Fring('^) to the function 

oo „ 

Fring{v) = J2<1 DiX)z{C)I,,^,{C) (67) 

9=1 



21 



is bounded by 



i^n„g(^) dk{pVik)ff^'\u)f^'\u) (68) 

where the inequahty ()25p has been also used. In view of the Lemma 1, it 
is clear that f^'^\y) (jMjl is finite for u < and thus -Fi.ing(z/) is also well 
defined in this extended range of chemical potentials. 

The ring contribution is expected to be small for a scaled potential 
when 7 — *■ 0. Indeed if the Fourier transform V^7(k) = K(k/7) of is 
introduced in (jfiSj) and the k integration variable is changed into k = 7p we 
find immediately 

Fring,,(^) < y / rfp(/3\>(p))Vf (^) = (^(7^), ^ < ^- (69) 

In (El, f\^Ki^) are still dependent on 7 through the single polymer partition 
function i^-y{q). 

If one introduces the scaled potential in the full expression (jU^ (setting 
k = 7p) one observes the subtle dependence of the parameter 7. There is an 
overall 7^ prefactor that manifests the smallness of the scaled potential as in 
the classical case. A 7 factor occurs also in the phase in the integrands of 
(I62|) and (j63j) in the combination 7A, representing the ratio of the thermal 
wave length and the potential range. The importance of these phases will 
therefore depend of the value of this ratio. Finally, 7 appears in the loop 
activity z^{C) in (jUHj) . For 7 small and 7A ^ 1 we can replace the above 
mentionned phase factors by 1 and approximate 

This will be the dominant term in the expansion of /cycle (^) since terms with 
A; > 1 cycles will obtain an overall 7'^'^ prefactor. 



6 Concluding remarks 

The use of quantum Mayer graphs proved useful in discovering the existence 
of an implicit equation ()49|) defining the density of an interacting Bose gas 
as function of its chemical potential. The knowledge of the precise form of 
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the equation requires the resummation of all multiply connected diagrams 
including the single vertex contribution. It should be recalled that in a 
series of papers, Lee and Yang have developped a diagrammatic formalism 
for the quantum statistical mechanical many-body problem that enables to 
calculate the thermodynamical quantities in terms of Boltzmann-type Ursell 
functions together with rules taking quantum statistics into account (see [21] 
for the general formalism and [22] for application to bosons). Working in 
the occupation number representation in momentum space, these authors 
also obtain a formally exact integral equation for the average occupation of 
modes, and in particular for the condensate density. The latter equation (a 
generalization of the Bogoliubov condensate equation) was shown to be exact 
in the thermodynamic limit by Ginibre [23], and it has also been established 
in the framework of infinitely extended states of Bose systems [24]. Our 
formulation is different in the sense that it provides a closed equation for the 
density of the interacting gas at fixed chemical potential that is particularly 
well adapted to the exploration of the neighborhood of the mean field limit. 
Moreover it provides an interesting link with the theory of polymers. 

Suppose that the resummation of multiply connected diagrams leading 
to ()49p can been performed. It is then not excluded that the resulting non- 
perturbative equation remains valid beyond the radius of convergence of the 
Mayer series. If so, the question of the Bose-Einstein condensation in an 
interacting gas could be examined on the basis of (|49|) . At least one example 
supports this hope: such an extension to the transition region agrees with 
rigorous results in the case of the mean field limit. This fact motivates the 
study of the self-consistent relation ()49|) beyond the mean field limit. We 
formulated this problem here using the scaling Q of the pair potential. In 
order to derive the form of function defining the relation (j49j) for the 
scaled potential one needs to know the small 7 asymptotics of 

• (i) the mean value of the Boltzmann factor of a single polymer (or 
normalized partition function) 

• (ii) the sum of ring diagrams representing the contribution of interact- 
ing polymers 

In the study of point (i) we could not use directly the known results (jHH) 
because of equal time condition imposed by quantum mechanics on interac- 
tions between different elements of polymers (see ((201) )• The existence and 
localization of the Bose-Einstein condensation depends in a crucial way on 
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the behavior of (i) for extended polymers. This however remains an open 
problem. Our paper provides only exact upper and lower bounds which, 
what is interesting, turn out to be qualitatively compatible with the classical 
result (jnU)- There are some indications that in the 7 — > limit the quantum 
calculation approaches the classical one. However, this question is at present 
not yet understood. We based the analysis of point (ii) on the remarkable 
fact that quantum Mayer ring diagrams can be summed in closed analytical 
form. We thus arrived at an analytic expression ()62|) for the sum of the rele- 
vant ring diagrams. However here again the small 7 limit of ()62j) involves the 
normalized single polymer partition function whose asymptotics for extended 
polymers remains to be derived. We hope that the continuation of the study 
of quantum Mayer graphs along these lines and progress in the understand- 
ing of polymer partition functions will eventually lead to the determination 
of functional relation beyond the mean field theory. 
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Appendices 

Appendix 1 : Proof of lemma 1 

The average energy (j37j) is calculated from its definition (introducing 
also the Fourier transform V(k) of the potential) 

{U), = Aq)-q^ (71) 

where 

J{q) = Dg{X) dsi ds2~S{s, - S2)V{X{X{s,) - Xis^))) 

= dkV{k) J' dsi J' ds26{si - S2) J Z}g(X)e'^'^[^(^i-^(^2)] 

(72) 

^From the basic rules for Fourier transforms of Gaussian measures we have 
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using the covariance (jl8|) 

J Dg(X)e'^''-[^(^^~^(^2)l = exp 



where 



(73) 



(74) 



Since both S{s) and Cq{s) can be considered as periodic functions of period 
q the double time integral in (f7^ reduces to 



dsi I ds2S{si — S2) exp 



2 ' 



q-l 

exp 

n=0 



-Cq{n) 



q-l 

exp 

n=l 



Cq{n) 



q (75) 



When this is introduced in (|7^ and (ffTjl (noting that / (iky(k) = \^(0)) one 
obtains finally 



{U)q = q dkV{k)g{k,q) 



with 



^(^. = 9 II exp 



2l2 



X'k 



-Cq{n) 



q/2 

5] exp 

71=1 



2u2 



X'k 



■CJn) 



(76) 



(77) 



The last equality results from the symmetry Cq{n) = Cq{q — n). If n < q/2, 
one has obviously Cq{n) > n/2 which leads to the bound 



q/2 

9{k, q) <Y1 

n=l 



2u2 



X'k 



-n 



< 



exp(^) - 1 



hence 



{U)g < qu+, u+ 



\/(k) 



exp(^; 



(78) 



(79) 



thus proving Lemma 1. 

It is instructive to exhibit the behaviour of in the scaling limit of a 
Kac potential V^7(k) = V^(k/7). Setting k = 7p in (f7^ gives 



exp(2^; 



-1 A2 
25 



4l'f,pM, -,^0 (80) 



Appendix 2 : Proof of lemma 2 

We denote {F)w = i^^(F|X(0) = 0) the normalized Wiener expectation 
of a functional -F(X) of paths X(t) starting from the origin at time t = and 
i?(F|X(ti) = Ri,X(t2) = R-2) the conditional Wiener expectation for paths 
starting in Ri at time ti and ending in R2 at time ^2- One has in particular 

{F)w = J dRE{F\X{0) = 0, X(t) = R), t>0 



exp(- 



|Ri~R2p 



E(l|X(ti) = Ri, X(t2) = R2) = (2^(^2-5')^/^ ' ^ ^^^^ 

The fact that Brownian motion is an homogeneous process implies the follow- 
ing symmetry relations under time and space translation, and space inversion: 

E(F|X(ti) = Ri,X(t2) = R2) = S(F,|X(ti + r) = Ri,X(t2 + r) =R2) 

= E(FR|X(ti) = Ri + R,X(t2) = R2 + R) 
= E(F„|X(ti) = -Ri,X(t2) = -R2) (82) 

with F.(X(-)) = F(X(- + r)), Fk(X(-)) = F(X(-) + R) and F_(X(.)) = 
F(— X(-)). Consider now the functional 



^gi,g2(X) = exp 



^ r ds r dt~6{s - t)(l - 5f^'j,[,^])V(X(.) - X(t)) 

Z Jqi Jqi 



(83) 

where the index (gi, ^2) means that the functional depends on the path only 
when t is in the interval [(3'i,(3'2]- In the present notation the normalized 
Brownian bridge average of a functional F reads 

J D^iX)F{X) = {2nqf/'E{F\X{0) = 0, X(g) = 0) (84) 

and therefore the vertex function K,{q) ()54|) for a q-particle loop X is 

K{q) = (27rg)3/2E(Fo,,|X(0) = 0, X(g) = 0) (85) 

In Fo giX) we suppress the interaction between the sets of particles 1, . . . , gi 
and gi + 1, . . . , g. Since V is positive this leads to the inequality 

Fo,,(X)<Fo,,,(X)F,,,,(X) (86) 
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implying in (|85p in view of the Markov property of the Brownian motion 
K{q) < {27rqf/^ JdRE{Fo,^,\X{0) = 0,X{q,) = R)E{F,,jX{q,) = K,X{q) 
= {27^qf'^ j t^RE(Fo,,jX(0) = 0,X(gi) = R)E(Fo,,, |X(0) = O.Xfe) 

with q = qi + q2- The second hne follows from the symmetry relations ()82p 
and the fact that the potential is invariant under space translations and 
inversion. Considering now the Wiener expectation of F^ qiX) we establish 
in the same way 

,gi+<?2)vK ^ {Po,qi)w{Po,q2)w (88) 

again as a consequence of the Markov property of the Brownian process and 
the invariance (jH2I)- We exploit the inequalities (jH7|) and (jSH|) as follows. We 
first use (jHTjl to relax the constraint of closed path, noting from ()81|) that 
£;(Fo,,jX(0) = 0,X(g2) = R) < (27rg2)~^/^ and from ^ 

/t(g) < ( 1 + ^ j (i"o,,i)w^, g = gi + g2, g2 > 1 (89) 

Then we set q\ = nqo, q2 = n, q = n{qo + 1) for some fixed integer go > 1- 
From (jH^ and the iteration of (jHHj) 

nin{qo + 1)) < (1 + qof'i{Fo,q,)wr (90) 
If one sets ^ 

= + 1) ^"^(^o^Jw > 0, r = go + 1 (91) 

one obtains the result of the lemma 2. 

Appendix 3: The Bose gas in an external potential 

Motivated by the experimental situation of Bosonic atoms in traps, we 
briefly show in this Appendix how the formalism works in presence of an 
external field (more details can be found in [25]). If a one-body external 
potential V^^^lr) is introduced, the formula for the grand canonical partition 
function S is still given by (jl4|) with the only change that the loop activity 
z{C) is replaced by 

~zic) = z(/:)e-^^^"(^) 

^ext(^) = rrfsr^"*(R + AX(s)) (92) 
Jo 
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Since the the density is non uniform, it is appropriate to consider here the 
average total particle number (N) 

(iV) ^ r'|lnH 

CO „ / n \ n 

= E ~ / ^-^1 ■ ■ ■ ^-^^ E n z{Ck)u{Ci, ...,Cn) 

n=l \£=1 / fc=l 

OO „ n 

= Et — dCi---dCnqiT\z{Ck)u{Ci,...,Cn) (93) 

^1 in -1)1 J ^Jl 

where we have introduced the Mayer expansion of InH and u{Ci, . . . , £„) is 
the Ursell function (jHH). 

The external potential will be chosen positive and confining with suffi- 
ciently fast growth at infinity so that / dr exp{—i3V^^^{r)) < oo. In order to 
allow for a well defined infinite particle number limit it is necessary to scale 
the external potential as 

\/J-t(r) = ^^"^(cjr) (94) 

with V^^^lr) a fixed positive confining potential. As the scaling parameter lj 
tends to 0, V^^*(r) becomes less confining and the average particle number 
{N){j3, n,u}) diverges. The proper quantity that remains finite in the limit is 
the product lj^{N){P, fi,uj). The experimental situation for atoms in traps 
appears to be well described in this asymptotic regime [1]. We have indeed 

Proposition 5 

In the range of convergence of the Mayer series for the uniform Bose gas, 
\\muj^QUj'^{N){l3, H^uj) = A/'(/3,/i) exists and 



X{P,fi) = j dvp{P,ii-V'-\v)) (95) 
with p{f3, yu) the density of the uniform Bose gas. 

In the scaling limit defined by (j94p . the gas is locally uniform at point r with 
a space dependent chemical potential /i(r) = /i — V^^^{r). 

We write V^^'^C) = ^^^^(cuR, q, cuX) and u{Ci, ...,£„)= u(Ri, gi, Xi, . . . , R„, g 
The spatial part of the integral in the n**^ order term of the series (j93j) reads 
(with the additional factor u^) 



a;3 



/n 
dK,dR2 • • ■ rfRn n e-'^^^"(-^-^--^'=)M(Ri, gi, Xi, R2, q^, X2, . . . , R„, X 
; — 1 



fc=l 
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k=2 

X M(0,gi,Xi,R2,g2,X2, . . . ,R„,gn,X„) (96) 

where one has made the sucessive change of variables R^ — > Rjt + Ri, k = 
2, . . . , n, then cjRi = r and used that the Ursell function is translation invari- 
ant with respect to its spatial variables. Since limj^^o V^'^^ir+uj'R.k, Qk, ^^k) = 
qi.V^^^{r) and the Ursell function is jointly integrable on R2, . . . , R^ the ex- 
pression tends by dominated convergence to 

Jdr f[ e-^'=^^"('-) [1 rfR2 ■ ■ ■ dKMO, Xi, gi, R2, X2, gs, • • • , Rn, X„, q„ 

(97) 

The first bracket, when combined with the loop activity z{C) ()22j). simply 
yields the shifted local chemical potential /i(r) = fi — V"'^^'(r) whereas the 
second bracket is the spatial part of the integral of the Ursell function oc- 
curing in the loop density series of the uniform gas evaluated at point 
Ri = 0. When this is introduced in the series for uj^{N){P, fi,uj) and 
the integrations on the internal degrees of freedom of the loops are taken into 
account (using also ((221) ), one obtains the result in the limit — > 0. 

To show convergence, one introduces the expresssion ()9(i|l in the series 
dnni) and majorize factors e-/3^'''*(r+'^R-fe.9fe,'^Xfe) A; = 2, . . . , n : 



n=l qi=l '' 

n „ 

X W z{Ck) / dC2 ■ ■ ■ dCnU{Cu £2, ■ • ■ , Cn) 

k=l 

00 ^ ^ 

= Y.<li D,,{X,) / rfre-'^^""('-'''^'-^^)pioop(^i) (98) 

where pioop(>Ci) is the Mayer series for the uniform system evaluted at 
£1 = (0, gi, Xi). Using again the positivity of V^^^^ijc) and Jensen inequality 
we have 

< j /^,,(Xi) j drexp(-/5^'rfsV;xt(r + c^AX(s)) 

< j Dg^{Xi) ds J drexp{-f3Ve^t{r + ujXX{s)) = J rfre-^^"'('")(99) 
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The convergence of Mayer series (j29p is established in Proposition 2 with 
bounds that are independant of the shape Xi of the loop L\. Combining 
this previous analysis with (j^ . one gets eventually 

uf'{N){^,y.,uj) < p(/i) / c/re-^^''^"-) (100) 



with p(/i) defined in Proposition 2. Since estimates are uniform with respect 
to uj, the existence of the limit (j^J^ follows again by dominated convergence. 

Moreover if the two-body potential is scaled according to 0, one obtains 
immediately from Proposition 3 the mean field limit for the trapped Bose 
gas 

Ar^f(/?,/i) = limAr(/3,/i,7)= / cirp^f(/5,/^-^'''(r)) (101) 

This result is established by other methods in [9] in a sligtly different situa- 
tion: the external potential considered in [9] has support in a box of volume 

and is scaled according to the size of the box (namely u = in (|94p ). 

The thermodynamics of the mean field trapped gas can be studied in 
detail from (jl01|) . If /i < pc = apo,c (the critical chemical potential of the 
homogeneous mean field gas), then /i(r) = fi — ^^^^(r) < fic for all r, so 
condensation does not occur anywhere. U fi > fic there is a local condensate 
of density [/i(r)/a — Po,c] at all points r such that p(r) > namely in the 
region of space A = {r |l^'^^*(r) < p — p^}. Then the total amount of 
condensate is 

A4ond(/3, l^)= [ dr [p(r)/a - po,c] (102) 

J A 

Obviously A/'cond(/3, p) vanishes ai fi = fic so the critical temperature Tc{Af) 
as function of the total particle number is defined as the solution of (jl01|) 
when fi = fif. and A/" is fixed, namely 

ACf (/5c, KPc)) = AT, p(/3) = apoAP) (103) 

As an example, for an harmonic potential V^^^lr) = 6|rp/2, one finds a criti- 
cal behaviour of the condensate fraction A/'cond(/3, A/')/A/' of the form (keeping 
now the particle number Af = A/mf (/?, p) fixed) [25] 

^^A,^ff\ T^T. (104) 

This is to be contrasted with the behaviours of the same fraction for the 
homogeneous free or mean field gas (~ 1 — (j^-) ) and for the free gas in 
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an harmonic potential (~ l — ^ijr^ — j ) when T approaches the corresponding 
critical temperatures. 

Clearly, when corrections to the mean field density of the homogeneous 
gas are known, they can be implemented in equation ()95|) and thus will also 
give interesting informations on imperfect trapped gases. 
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